Unified picture of Q-balls and boson stars via catastrophe theory 
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We make an analysis of Q-balls and boson stars using catastrophe theory, as an extension of the 
previous work on Q-balls in flat spacetime. We adopt the potential V3 (4>) = ^-<^ 2 — + A0 4 for 
Q-balls and that with [i = for boson stars. For solutions with \g rr — 1| ~ 1 at its peak, stability 
of Q-balls has been lost regardless of the potential parameters. As a result, phase relations, such as 
a Q-ball charge versus a total Hamiltonian energy, approach those of boson stars, which tell us an 
unified picture of Q-balls and boson stars. 
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I. INTRODUCTION 

Among non-topological solitons which appear in U(l)- 
symmetric scalar fields, objects existing even in flat 
spacetime are called Q-balls [ll-flTj]. while objects sup- 
ported by strong gravity are called boson stars [r3 - [l7| . 

Q-balls are typically supposed to be microscopic ob- 
jects. It has been argued that Q-balls generally exist in 
all supersymmetric extensions of the Standard Model and 
could play important roles in cosmology Q. For example, 
Q-balls can be produced efficiently in the Affleck-Dine 
mechanism and could be responsible for baryon asymme- 
try [i[ and dark matter Q . As for the stability of Q-ball 
in flat spacetime, analytic results were obtained under 
the thin-wall approximation |2| and under the thick-wall 
approximation [6[. A general criterion of stability was 
derived analytically by Pacceti Correia and Schmidt Q: 
a Q-ball is stable if and only if 



U — < 

Q duj 



(1.1) 



where ui and Q are the angular velocity of phase rotation 
and the Q-ball charge, respectively. In general cases, nu- 
merical calculation is necessary for having equilibrium 
solutions, and catastrophe theory [l8[ is a useful tool for 
finding their stability. Kusmatsev made a general argu- 
ment that catastrophe theory can be applied to the inves- 
tigation of molecules and solitons in various systems and 
found the criterion (jl.ip Q . In @ catastrophe theory was 
applied to two typical models of Q-balls to explore the 
stability for the whole parameter space, which includes 
the intermediate region between the thin-wall limit and 
thick-wall limit. In this paper, as well as in most of the 
previous works, stability means local stability, that is, 
stability against small perturbations. Absolute stability 
of Q-balls has been argued in [l(| ■ 
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On the other hand, boson stars have been studied as 
astronomical objects which could also contribute to dark 
matter [l2j]. For example, axidilaton star of ~ O.6M 
could account for a part of massive compact halo ob- 
jects [la ]. Supermassive boson stars of 10 6 -10 9 M Q have 
been discussed as an alternative to a black hole in the 
galaxy center [l4| ■ If we consider the evolution of boson 
stars in scalar-tensor theories, they could have gravita- 
tional memory: the strength of the gravitational constant 
at formation time could still be effective [l5j . As for the 
stability analysis of boson stars, catastrophe theory has 
also been used [la. 



We may note, in passing, that catastrophe theory has 
also been applied to the stability analysis of non-Abelian 
black holes 



Although the difference in theory between Q-balls and 
boson stars is solely the model parameters, the investiga- 
tions of their properties have been carried out separately. 
This is because Q-balls and boson stars have been dis- 
cussed in different contexts of particle physics or astro- 
physics. The purpose of the present paper is to obtain an 
unified picture of equilibrium solutions and their stabil- 
ity of Q-balls and boson stars. Our interest is essentially 
mathematical. However, gravitating Q-balls, or Q-stars, 
which are intermediate objects between Q-balls in flat 
spacetime and boson stars, have also been discussed [20I — 
[22]]. Therefore, an unified analysis of Q-balls and boson 
stars is also important for the study of astrophysical mod- 
els. 



This paper is organized as follows. In Sec. II, we derive 
equilibrium field equations and explain how to analyze 
stability using catastrophe theory. In Sec. Ill, we show 
numerical results of equilibrium Q-balls and boson stars 
and discuss their stability. In Sec. IV, we devote to 
concluding remarks. 
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II. ANALYSIS METHOD OF EQUILIBRIUM 
Q-BALLS AND BOSON STARS 

A. Equilibrium field equations 

We begin with the action 

S : / cf 4 x(£ G +£ ), 



£g = 



-9 



16ttG 

\sr»„. 



- V^g { —fTd^ ■ d u <p - v(<j>) } , (2.1) 



where <f> — (<f>i, (fa) is a SO(2)-symmctric scalar field 
and <j> = ' 4> = \/$L + 4%- We assume a spherically 
symmetric and static spacetimc, 

ds 2 = -a 2 (r)dt 2 + A 2 {r)dr 2 + r 2 (d6 2 +sin 2 6d(p 2 ). (2.2) 

For the scalar field, we assume that it has a spherically 
symmetric and stationary form, 



B. Models of Q-balls and boson stars 

As for Q-balls, which are present even in flat spacetime, 
we suppose the potential, 

9 

V 3 (<t>) := —<p 2 -n(f> 3 + \(t> 4 with to 2 , n, A > , (2.9) 

which we call V3 Model. Rescaling the quantities as [23j j 

to [a m 2 G 

t = mt, r = mr, ui = — , \i=—= — , K = — - — , 
m VAto A 



- 

= 9, 

m 

Q = XQ, 



(2.10) 



the field equations (|2.4j) . (|2.5j) and (|2.6j) with the poten- 
tial (12.91) are rewritten as 



(0i , ^2) = 0(r)(cos wt, sinwt). 
Then the field equations become 



^ ee A' + £(A 9 -1) 



(2.3) 



AirGrA 3 ^- 



i/2 2 
CJ 



2A 2 2c 



(2.4) 



ra 



a 



a' + g(l-A 2 ) 



-□c 



V 2 


W 2 2 


2A 2 


+ 2a 2 


a' 




Q 





y , (2.5) 



ujA 
a 



= A 2 ^. 



(2.6) 



where ' = d/dr. To obtain Q-ball solutions in curved 
spacetime, we should solve (I2.4j) - (|2.6[) with boundary 
conditions, 

A(0) = A(oo) = a(oo) = 1, 

A'(0) = a'(0) = <f/(0) = (f>(oc) = 0. (2.7) 

We also restrict our solutions to monotonically decreasing 
4>(r). Due to the symmetry, there is a conserved charge 
called Q-ball charge, 



d^xyf^gg^^idvfo - (fadvfa) = wl, 



Ar 2 6 2 

where 7 = 4tt / —dr. (2.8) 



(2.12) 



2 a' A 1 



jA 



,dV 



' + ( — ) ~4> = A 2 ^—. (2.13) 
r a A ] \ a 



For reference, we recall the parameter regions of w 2 
where solutions exist in flat spacetime. In this case, the 
field equation is 



d: 



2d0 -27, W 

-T3T = ---T7-W - 

dr l r dr dd 



(2.14) 



This is equivalent to the field equation for a single static 
scalar field with the potential = V — tu 2 (j) 2 /2. Equilib- 
rium solutions satisfying boundary conditions (|2.7I) exist 
if min(VL) < V(0) and d 2 K,/d0 2 (O) > 0, which is equiv- 
alent to 



i-y<* 2 <i- 



(2.15) 



The two limits uj 2 — > 1 — 4- and uj 2 — > 1 correspond to the 
thin- wall limit and the thick- wall limit, respectively. As 
we shall explain below, qualitative features of solutions 
change at ft, 2 = 2. 

On the other hand, as for boson stars, we assume the 
potential (|2.9[) with fi = 0. Their equilibrium solutions 
do not exist in flat spacetime. 

The condition (j2 . 15[) cannot apply to gravitating solu- 
tions. Actually, boson stars exist. However, as we shall 
show below, the criterion w 2 < 1 is valid even for gravi- 
tating Q-balls and boson stars. 
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C. Stability analysis method via catastrophe 
theory 

Let us discuss how we apply catastrophe theory to the 
present Q-ball or boson star system. An essential point is 
to choose behavior variable(s), control parameter's) and 
a potential in the Q-ball or boson star system appropri- 
ately. 

In [9j we argued that the total energy of the scalar 
field, 



d 6 x 



2m2 



k'\2 



V 



(2.16) 



is appropriate for a potential because the variation of E$ 
under fixed Q, 5E^/84>\q = 0, reproduces the equilibrium 
field equation (|2.14p . This is on the analogy of a poten- 
tial in a mechanical system, where the potential F{x) is 
defined in such a way that dF/dx — at equilibrium 
points. 

A nontrivial issue in curved spacetime is the choice 
of the corresponding total energy since there are many 
definitions for total energy. However, we can conclude 
that the Hamiltonian energy E, which we shall calculate 
below, is the appropriate because 8E /54>\q = SE/Sg^ = 
0, reproduces the equilibrium field equations (I2.4j) - ([2.6p . 

From the Lgrangian (|2.ip with the coordinate system 
(|2.2p , one finds the canonical momentum of <f> and the 
matter part of the Hamiltonian density, 



V 



1 — g = aAr sin ( 



(2.17) 



3i<f> = 3'a(f>a—£4> — 




(2.18) 



where ' = d/dt. Under the stationary condition 
we obtain 



'L0 2 <f> 2 (0') 2 t . 



2a 2 



2A 2 



-gp<p- (2.19) 



Similarly, we can consider the canonical momentum of 
the spatial metric hij and the gravity part of the Hamil- 
tonian density, 



dh t 



(2.20) 



3iV 



7T lJ hij - £g 



9. 



16ttG 




(2.21) 



where we have used the static condition, hij — 0. Using 
one of the field equations, G\ = SnGp,/,, we obtain the 
total Hamiltonian, 



E : 



d 3 x(% G + J{ ) = lim 



rV 
2GA' 



(2.22) 



If we define the gravitational mass M by the asymptotic 
behavior of the metric, 



Eq. ([2~22| reduces to 



2GM 



as r — > oo, 



2 



(2.23) 



(2.24) 



In the previous work [20( the total Hamiltonian was cal- 
culated as E = M; however, what we have shown is that 
the correct formula is (|2 . 24[) . We also use the normalized 
quantity 



E 



X 



-E. 



(2.25) 



Because the charge Q and the model parameter (s) of 
V{4>) can be given by hand, they should be regarded 
as control parameters. In flat spacetime, V3 Model es- 
sentially has only one parameter, p? . In curved space- 
time, on the other hand, the normalized gravitational 
constant k is another control parameter, which represents 
the strength of gravity. 

To discuss a behavior variable we consider an one- 
parameter family of perturbed field configurations 4> x (r) 
near the equilibrium solution <p(r). Because dE[4> x ]/dx = 
(SE I '5(f>x)d(j) x I dx = when <p x is an equilibrium solution, 
x is a behavior variable. Although an explicit choice for 
x is not unique, we choose uj 2 as a behavior variable. 

According to Thorn's theorem, if the system has two 
control parameters in a mechanical system, there is essen- 
tially one behavior variable; if the system has three con- 
trol parameters, there are one or two behavior variables. 
Because the present Q-ball or boson star system con- 
tains (Q,/i 2 ,K), we speculate that each has two behavior 
variables, oj 2 and 4>(0), and falls into hyperbolic umbilical 
catastrophe. However, because the stability structure of 
equilibrium solutions in three-parameter space (Q, /I 2 , k) 
is very complicated and our interest is how gravitational 
effects change the stability structure, in the following, we 
discuss the stability structure of equilibrium solutions in 
two-parameter space (Q,n) under fixed jl 2 . 

Our method of analyzing the stability of Q-balls and 
boson stars is as follows. 

• Fix the value of p, 2 . 

• Solve the field equations (|2.4|) - (12.6[) with the bound- 
ary condition (|2.7p numerically to obtain equilib- 
rium solutions 4>(r) for various values of uj and k. 

• Calculate Q for each solution to obtain the equilib- 
rium space M = {(x, Q, k)}. We denote the equa- 
tion that determines M by f(x,Q,K) = 0. 

• Find folding points where dQ/dx = or dn/dx = 
0, in M, which are identical to the stability-change 
points, E = {(x, Q, k) \ df/dx = 0, / = 0}. 
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FIG. 1: (a) Q-E relation and (b) Q-Cj 2 relation for jl 2 = | in 
Vz Model. We compare solutions for k = 0, 0.006 and 0.6. In 
the case of gravitating Q-balls (k 7= 0) cusp and spiral struc- 
tures can be seen in (a) and (b), respectively. The maximum 
of Q and the local minimum value of Q for k = 0.006 are 
labeled as A and B, respectively. We interpret solutions with 
solid (dotted) lines stable (unstable). 

• Calculate the energy E by (I2.22[) for equilibrium so- 
lutions around a certain point in E to find whether 
the point is a local maximum or a local minimum. 
Then we find the stability structure for the whole 
M. 



III. STABILITY OF GRAVITATING Q-BALLS 
AND BOSON STARS 

Let us consider Q-balls of V3 Model (12.91) . In flat space- 
time (k — 0), stability structure falls into two classes, 
jl 2 < 2 and fi 2 > 2 @: 

• jl 2 < 2: All equilibrium solutions are stable. 



• jl 2 > 2 : For each /2 2 , there is a maximum charge, 
Q m ax, above which equilibrium solutions do not ex- 
ist. For Q < Qmax, stable and unstable solutions 
coexist. 

A. Gravitating Q-balls for fj, 2 < 2 

In this subsection, we fix p 2 = | as an example of 

jl 2 < 2. Figure [T] shows a plot of Q versus E and that of 
Q versus ui 2 for equilibrium Q-ball solutions. In the case 
of k = there is one-to-one correspondence between Q 
and E while cusp structures appear in the case of k =/= 0, 
as shown in (a). The maximum of Q (labeled as A for 
k = 0.006) and the local minimum (labeled as B for 
k = 0.006) appear in the case of gravitating Q-balls. At 
the point £?, another cusp strcture appears and is far 
smaller than that at the point A. This sequences of cusp 
structure continue and we stopped calculation where the 
4th cusp structure appears. Similar structures have been 
reported in [13, [H|. The Q-maximum for n = 0.6 is 
far smaller than that for k = 0.006. Figure QJb) shows 
that the stability criterion (jl.ll) . which was established 
by perturbation on equilibrium Q-balls in flat spacetime, 
cannot apply to gravitating Q-balls. 




FIG. 2: <f> as a function of r for k — (dot-dashed lines) and 
for ft = 0.006 (solid lines and a dotted line). For each k, we 
take Q = 200 and 800. 

To see how gravity changes stability structure of equi- 
librium solutions, we show the profiles of the scalar field 
for k — and those for n = 0.006 in Fig. [2] We find 
that gravitational effects reduce the Q-ball size and they 
become larger as Q becomes larger. We also find that 
there are two solutions for Q — 800 and k = 0.006. The 
solution with larger cf> near the origin corresponds to the 
dotted line in Fig. [TJ We suppose that Q-balls with dot- 
ted lines cannot support itself and will collapse or dis- 
perse and the point A would be the point where stability 
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transformation uj — > — uj changes nothing but Q — > — Q. 
Stability of the solutions is determined by their energy 
E, as calculated in Fig. [TJ 




FIG. 3: Structures of the equilibrium spaces, M = 
{(w 2 , k, Q)}, and their catastrophe map, x(M), into the con- 
trol planes, C = {(jl 2 ,Q)}, for jl 2 = |. Blue lines and red 
lines in M represent stable and unstable solutions, respec- 
tively. In the regions denoted by S, SU and N on C, there 
are one stable solution, one stable solution and one or more 
unstable solutions, and no equilibrium solution, respectively, 
for fixed (k, Q). 



changes. Accordingly, there also exists the maximum of 
the Q-ball size due to gravity, as was pointed in [22| . We 
thus reasonably understand that the maximum of Q for 
K = 0.6 is far smaller than that for k — 0.006. 

Figure [3] shows the structures of the equilibrium spaces, 
M = {{Cj 2 , k, Q)}, and their catastrophe map, x(M), into 
the control planes, C — {{p?,Q)}, for ft 2 = |. x(M) 
shows that in the regions denoted by S, SU and N on C , 
there are one stable solution, one stable solution and one 
or more unstable solutions, and no equilibrium solution, 
respectively, for fixed (n,Q). For example, for k = 0.6, 
which is the case chosen in Fig. [TJ we can confirm that 
unly a stable solution exists below Q ~ 0.5 while one 
stable solution and one or more unstable solutions exist 
in the region from Q~0.5toQ~0.9. Here we demon- 
strate only the results for w > 0, that is, Q > 0. This 
does not mean Q and Q are always positive; the sign 




FIG. 4: Stability interpretation via catastrophe theory for 



P 2 



| for the gravitating case with k = 0.006 (The qualita- 



tive properties are same for other k). The lines when we fix 
the control parameter Q in (a) are depicted by the quadratic 
curves in (b). ai and 62 are regarded as the potential mini- 
mum and the maximum in (c), respectively. 

We explain our interpretation via catastrophe theory 
by exemplifying the case with jj? = | and n = 0.006 
(Qualitative properties are not changed for other k.). We 
identify Q =const. lines in Fig. [4] (a) with the quadratic 
curves in Fig. U (b). If we adopt the view point that 
stability changes at the point A as we mentioned above 
and observe Fig. 0] (b), we notice that (—^(0)) is more 
appropriate for a behavior variable than uj 2 . Then, as 
we show in (c), 02, 62 and A can be interpreted as the 
potential minimum, maximum and the inflection point, 
respectively. 

This case can be understood using the fold catastrophe 
f(u) = u 3 + tu where u and t are the bahavior variable 
identified with (—0(0)) and the control parameter iden- 
tified with Q, respectively. 

We should reveal what causes the difference from the 
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FIG. 6: Q-E relation for ft 2 = 2 in V 3 Model. We compare 
solutions for k = and 0.05. In the case of gravitating Q- 
balls [k 7^ 0) cusp structures can be seen as in the case with 
ti 2 - 5 



FIG. 7: (a) Q-E relation and (b) Q-Cj 2 relation for fi 2 = 5. 
Intrinsic difference can be seen for = 0.3. 



flat case. We show behavior of the metric A for the so- 
lutions with Cj 2 = 0.098, 0.16, 0.27, 0.63 and 0.98 for 
k = 0.006 in Fig. [5] We should notice its relation to Fig.0] 
(b). Naively speaking, solutions having larger \ A— 1| at 
its peak have larger (—0(0)). We have confirmed that, 
independent of k, solutions having \A— 1| ~ 1 correspond 
to solutions expressed by dotted lines in Fig. [TJ There- 
fore, we can suppose that the intrinsic difference from the 
flat case can be characterized by \A— 1|. 

We have also confirmed that these properties can be 
seen for other fi 2 . As an example, we exhibit Q-E re- 
lation for fi 2 — 2 in Fig. 6. In the case of gravitating 
Q-balls (k ^ 0) cusp structures can be seen as in the 
case with fx 2 = |. Other diagrams and related proper- 



B. Gravitating Q-balls for fi 2 > 2 

Next, we fix fi 2 — 5 as an example of jl 2 > 2. Figure [7| 
shows a plot of Q versus E and that of Q versus w 2 for 
k = 0, 0.2 and 0.3. In flat spacetime (k = 0), contrary to 
the case of jj 2 < 2, Q has a maximum Q max , and for each 
Q{< Qmax) there is one stable solution and one unstable 
solution The solution of the Q-maximum coincides 
with that of the ^-maximum. Figure [7] indicates that 
gravitational effects for the case of k = 0.2 is not so large, 
compared with the case of n = 0.006 and ft 2 — | in Fig. 
[TJ This is simply due to the smallness of the gravitational 
mass IE. 

In the case of k = 0.3, we find a cusp structure sim- 
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FIG. 8: Structures of the equilibrium spaces, M= {(cl> 2 , K, Q)}, 
and their catastrophe map, x(M), into the control planes, C = 
{(ft 2 , Q)}, for fi 2 — 5. Blue lines and red lines in M represent 
stable and unstable solutions, respectively. In the regions 
denoted by S, SU and N on C, there are one stable solution, 
one stable solution and one or more unstable solutions, and 
no equilibrium solution, respectively, for fixed (k,Q). 



ilar to that for p? = |. We call it a degenerate cusp 
structure. We also confirmed that solutions written by 
a dotted line for k — 0.3 have the metric \A — 1| ~ 1 at 
its peak. Then, again, we suppose that solutions with 
strong gravity \A — 1| ~ 1 show similar cusp structure 
regardless of the potential parameters. 

In addition, we have solutions with another sequences 
of cusp structures at far larger Q and E around Q ~ 3.5. 
We call the former sequence the low energy branch, while 
the latter sequence the high energy branch. The low en- 
ergy branch is similar to the sequence of p? = | in Fig. 
[TJ which suggests that the stability structure is also sim- 
ilar. On the other hand, we suppose that any solution 
in the high energy branch is unstable since it basically 
corresponds to that written by a dotted line for k = 0, 
which is already shown as unstable. In the regions de- 
noted by S, SU and N on C, there are one stable solution, 
one stable solution and one or more unstable solutions, 
and no equilibrium solution, respectively, for fixed (k, Q). 
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FIG. 9: Stability interpretation via catastrophe theory for 
fi 2 — 5 for the flat case. The lines when we fix the control 
parameter Q in (a) are depicted by the quadratic curves in 
(b). The points 02 (02) are supposed to be stable (unstable) 
since it is expressed as the potential minimum (maximum) 
shown in (c). 



Figure [5] shows the structures of the equilibrium spaces, 
M = {(tl> 2 , k, Q)}, and their catastrophe map, x(^-)> m t° 
the control planes, C = {(p?, Q)}, for jl 2 = 5. 

We see that one sequence diverges into two branches at 
k w 0.28, above which each branch exposes spiral struc- 
tures. The upper branch in Fig. [5] is analogous to the 
solution sequence for jl 2 = -| in Fig. [3] This also indi- 
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cates that, as the strength of gravity becomes large, the 
dependence on the potential shape diminishes. 
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FIG. 11: Q-E relation for ft 2 = | which is qualitatively same 



(K0) 



as the case with /i 



near the point A. As for the high energy branch, since 
the point C2 corresponds to that in Fig. El (c), we can 
naturally suppose that this branch is unstable. 

For completeness, in Fig. [Til we a l so show Q-E relation 
for jl 2 = | with k = and 0.0792 which is qualitatively 
similar to the case with jl 2 = 5. The solid and the dotted 
lines correspond to the stable and the unstable solutions, 
respectively. As for the case ft, 2 = I with n — 0.0792 
the only difference from the case jl 2 — 5 with k = 0.3 is 
that there are stable solutions in the high energy branch. 
However, this is not the intrinsic difference. Actually, for 
the case jl 2 = 5 with k ~ 0.28, there are stable solutions 
in the high energy branch which correspond to the point 
a2 as in the low energy branch. As we can expect from 
this diagram, other relations also resemble to those with 
A 2 



-■ 2 = 5. 



FIG. 10: Stability interpretation via catastrophe theory for 
jl 2 = 5 for the gravitating case with k = 0.3 (The qualitative 
properties are same for other k > 0.28. The case for k < 0.28 
can be understood as in the same way in the flat case.). The 
lines when we fix the control parameter Q in (a) are depicted 
by the quadratic curves in (b). In this case, potential plane 
is interpreted as Fig. [4] (c) . 

We can understand its catastrophe type using fold 
catastrophe as the case with p? = |. We also iden- 
tify Q =const. laines in Fig. |H] (a) with the quadratic 
curves in Fig. [9] (b) where flat solutions are shown. In 
this case, there is a one-to-one correspondence between 
uj 2 and 0(0). Therefore, both variables can be behavior 
variables as shown in Fig. (c). 

We also exhibit the corresponding figures with k = 0.3 
in Fig. QJJJ uj 2 -(~ (f)(0)) relation changes drastically as 
shown in Fig. [10] (b). This figure suggests that we should 
use (—(f>(0)) as a behavior variable similar to Fig. 0] (c) 



C. Boson stars 

Now we discuss boson stars with the potential (|2.9[) 
with fi = 0. Figure [T2l shows plots of (a) Q-E, (b) Q-Cj 2 
for equilibrium solutions of boson stars. Degenerate cusp 
and spiral structures are seen as in the case of gravitating 
Q-balls for jl 2 = | and those for fi 2 = 5 with k > 0.28. 

Figure [T3l shows the structures of the equilibrium spaces 
and their catastrophe map x(-M-) into the control planes 
for jj? = 0. In the regions denoted by S, SU and N 
on C, there are one stable solution, one stable solution 
and one or more unstable solutions, and no equilibrium 
solution, respectively, for fixed (n,Q). We see that qual- 
itative characteristics of the equilibrium space and its 
catastrophe type are the same as those for fi 2 =5/3. We 
have confirmed \A — 1| ~ 1 at its peak in the solutions 
corresponding to the spiral curves or near the stability 
change points. We therefore conclude that, if gravity is 
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and boson stars. 

Therefore, if we discuss the possibility of Q-balls or 
boson stars as dark matter candidates, our work would 
be useful. This work should also be extended to the V4 
model which will appear in our companion paper. 



50 100 150 200 250 300 




FIG. 12: (a)Q-E, (b)Q-£l> 2 relations in the model of boson 
stars for k = 0.006. 



so strong as \A — 1| ~ 1 at its peak, catastrophic struc- 
tures of Q-balls approach those of boson stars, regardless 
of the potential shape. 



FIG. 13: Structures of the equilibrium spaces, M — 
{(ti 2 , k, Q)}, and their catastrophe map, x(M), into the con- 
trol planes, C = {(/i 2 , Q)}, for fl 2 = 0. In the regions denoted 
by S, SU and N on C, there are one stable solution, one stable 
solution and one or more unstable solutions, and no equilib- 
rium solution, respectively, for fixed (k, Q). 



IV. CONCLUSION AND DISCUSSION 

We have reanalyzed stability of gravitating Q-balls for 
a V3 model and boson stars for a V3 model with [i = 0. 
For solutions with \g rr — 1| ~ 1 at its peak, stability of Q- 
balls has been lost regardless of the potential parameters. 
As a result, phase relations, such as Q-E, approach those 
of boson stars, which tell us an unified picture of Q-balls 
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For this reason, present quantities are written using the 
previous quantities as 



M = — , 
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As for the total energy E, care should be taken. In our 
previous paper, only H$ has been considered while + 
Ho are considered, in the present paper, to define the 
total energy for including the gravitating case. Because of 
this difference, in flat case, the total energy in the present 
paper is half of that in the previous paper. Therefore, we 
obtain the relation as 
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